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1. Introduction 

In recent years the standards were established for the perturbative QCD descrip- 
tion of various characteristics of hadron jets produced in e + e~ annihilation. These 
standards include: 

• all-order resummation of double- (DL) and single-logarithmic (SL) contribu- 
tions due to soft and collinear gluon radiation effects, 

• two-loop analysis of the basic gluon radiation probability and 

• matching the resummed logarithmic expressions with the exact O (a^) results. 

Such programmes were carried out for a number of jet shape observables such as 
Thrust (T) and heavy jet mass (M H ) [1], C-parameter [2] and jet Broadenings (total 
B T and wide-jet broadening By/) [3]. 

Perturbative description of jets produced in processes other than e + e~ annihila- 
tion poses more difficulties, as the structure of jets (both average jet characteristics 
and distributions) becomes sensitive to details of the underlying hard interaction. 
For example, characteristics of quark jets produced in the current fragmentation re- 
gion in Deep Inelastic Scattering (DIS) depend on the Bjorken x. In spite of these 
complications a steady progress is being made in this domain as well [4]. 

Given a (relative) perfection of perturbative QCD technologies, it became possi- 
ble to aim at deviations of the measured hadronic characteristics from the correspond- 
ing perturbative predictions, in a search for genuine non-perturbative confinement 
effects. As is well known, these deviations, for a broad variety of jet observables, 
amount to sizable 1/Q corrections, with Q the hardness scale (the total e + e~ anni- 
hilation energy), for reviews see [5,6]. 

As far as physics of e + e~ annihilation is concerned, till now these developments, 
both in the perturbative and non-perturbative sectors, were confined to two-jet en- 
sembles which constitute the bulk of events. Little (if any) attention has been paid 
to multi-jet ensembles, in particular to three-jet events. In spite of being obviously 
more rare [o^jo t ot oc <y^~ 2 , with N the number of jets), selected multi-jet config- 
urations are of special interest as they, so to speak, are subject to more quantum 
mechanics than unrestricted (mainly two-jet) hadron production events. 

Indeed, in the latter case it is known that the gross inclusive features of particle 
production can be described in probabilistic terms by imposing angular ordering [7] 
on successive 1 — > 2 intra-jet parton decays. It suffices to implement strict an- 
gular ordering, 9 i+ i < 9i, the proper running coupling, a s (k±), and the standard 
parton splitting probabilities to ensure the next-to-leading accuracy of the perturba- 
tive description of inclusive energy spectra, mean multiplicities and the multiplicity 
distribution [8]. 
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Thus, a transparent and powerful probabilistic technology exists for predicting 
(with next-to-leading accuracy) inclusive observables. At the same time, the very 
selection of events (say, three-jet events) necessarily makes a measurement "less 
inclusive" and destroys, generally speaking, the probabilistic picture, and one's hold 
on the SL accuracy with it. 

Of special interest for two-jet configurations is the kinematical region of two 
narrow jets (1 — T, Mfj/Q 2 ,C, B T W <C 1) where multiple radiation (Sudakov sup- 
pression) effects are essential. An analogue of this region for three-jet events is 
near-to-planar kinematics. Three hard partons — qq and a gluon g — form a plane. 
Secondary gluon radiation off the three-prong QCD antenna brings in aplanarity. In 
what follows we choose the out-of-plane transverse momentum K ont as an aplanarity 
measure. The distribution of events in K out is subject to DL suppression in the region 
of small aplanarity, K out /Q <C 1, where normal accompanying radiation is vetoed. 

The corresponding suppression factor, and thus the resulting i^ out -distribution is 
easy to predict in a DL approximation which takes care of the leading contributions 
O {[a s ln 2 ((5/i^out)] n ) i n & U orders while disregarding subleading SL corrections of 
the order of [a s ln(Q/K ont )] m . The answer is given simply by the product of three 
proper QCD Sudakov factors, F 2 ■ F g , which veto radiation of gluons with out-of- 
plane momentum components exceeding K out , off the three hard partons treated as 
independent emitters. 

The DL approximation is known to be too rough to be practically reliable. Im- 
proving it proves to be nontrivial a quest: at the level of SL terms the geometry of 
the underlying three-jet configuration enters the game: in addition to intra-jet par- 
ton multiplication one has to take into consideration inter-jet particle production. 
The latter, however, does not admit "classical" probabilistic interpretation: gluon 
radiation in-between jets results from the coherent action of all three elements of the 
antenna. 

This is a general feature which complicates the analyses of other, more simple, 
characteristics of three-jet ensembles as well. For example, contribution of coherent 
inter-jet particle flows enters, at the level of the next-to-leading O (\fa^) correction, 
into perturbative prediction of the mean particle multiplicity in three-jet events, 
making it event-geometry-dependent [9]. 

In this paper we attempt, for the first time, the all-order perturbative analysis 
of the -Kout-distribution in three-jet e + e~ annihilation events, aiming at SL accu- 
racy. In what follows we will single out and resum logarithmically-enhanced DL 
and SL contributions and systematically neglect relative corrections of the order 
0(a s ). The latter belong to the non- logarithmic normalization factor ("coefficient 
function") 1 + ca s + . . . , whose first coefficient, c, can be found by comparing the ap- 
proximate resummed result with numerical calculation based on the exact a 2 matrix 
element [10]. 

To accommodate all essential SL contributions one has 
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1. to take into account soft inter-jet gluon radiation, 

2. to analyse corrections due to hard intra-jet parton decays, 

3. to define an event plane for a multi-parton system and properly treat kinemat- 
ical effects due to parton recoil, 

4. to prove, for three-jet environment, soft gluon exponentiation (at the two-loop 
level) and the prescription for the argument of the running coupling that enters 
the basic gluon emission probability. 

In the present paper we address these issues. 

The answer we derive has the following key features: 

• Kinematical constraints which determine an event plane are rather complicated 
but can be resolved with a help of multiple Fourier-Mellin representation which 
allows for exponentiation of multiple radiation in the parameter space. 

• Multiple soft radiation off the three-parton system can be resummed and ex- 
ponentiated in terms of three colour dipoles that together determine the colour 
structure of (and accompanying particle flows in) the event. 

• For the cumulative i^ out -distribution (as well as for other sufficiently inclusive 
observables) inclusive treatment of the two-parton decay of a gluon emitted by 
the qqg system results in the running of the coupling constant. 

• The running a s which describes the intensity of gluon emission off each dipole 
(qg, qg, qq) depends on the invariant transverse momentum of the gluon with 
respect to two partons that form the corresponding dipole. 

• The exponent can be represented as a sum of three basic parton "radiators" 
each of which describes one-gluon emission off a single hard parton, weighted by 
the colour factor of this parton. This radiation can be treated as independent 
provided a proper hardness scale is ascribed to the parton radiator. 

• Essential SL contribution ("hard" intra-jet and coherent inter-jet corrections) 
can be conveniently embodied into the scales. Effects of the large-angle soft 
(inter-jet) radiation make these scales event-geometry-dependent and different 
for each of the three primary partons. 

• The structure of the hardness scales of the parton radiators entering in the total 
fT out -distribution has a clear geometrical interpretation: the scale Q a for the 
parton a is proportional to the invariant transverse momentum of this parton, 
p ta , with respect to the hyper-plane formed by the other two hard partons. 
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• Equivalent expressions for the radiators of the total -K" out -distributions can be 
constructed which smoothly interpolate between 3- and 2-jet configurations. 

As has been already said, kinematics of three-jet observables complicates the analysis. 
As a result, the final expressions are rather cumbersome as they involve 4- and 5- 
dimensional integrals. 

The kinematical constraints, and thus the final formulae, are somewhat simpler 
for the distribution of K ont accumulated in the right hemisphere, i.e. the one that 
has the smallest transverse momentum with respect to the thrust axis (for a review 
see [11]). 

Quantitative analysis of the predictions and numerical results are discussed in a 
separate publication [10]. 

The paper is organised as follows. 
In Section 2 we define the aplanarity measure K OVLt , derive kinematical relations 
defining the event plane and discuss resummation of soft gluon radiation for near- 
to-planar 3-jet events. 

Section 3 is devoted to perturbative analysis of the -K" ut-distribution in the right 
hemisphere. Here we develop technique for analysing and embodying all necessary 
SL corrections into the all-order perturbative result. 

In Section 4 we apply this technique to derive, with SL accuracy, the perturbative 
prediction for the total -K" ut-distribution in 3-jet events with given kinematics (thrust 
T and thrust-major T M )- 

Technical details are confined to Appendices. 



2. Aplanarity and soft parton resummation 

We consider e + e~ annihilation events with almost planar configuration of final state 
particles. Such event are characterised by the inequality 

T~T M >T m , (2.1) 

with T the thrust, T M and T m the so-called thrust-major and thrust-minor. We 
study the distribution of T m . We shall refer to e + e~ events in this phase space region 
as 3-jet events. 

Thrust is defined as 

TQ = max <^ ^ I ™Ph I f = I ^rPh | = ^ [ Phz | , (2.2) 
I h J h h 

where the sum runs over all particles h produced in a given event with the total 
center of mass energy Q. Hereafter we choose the z-axis to lie along the thrust axis, 
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Or- The thrust-major is defined analogously; it maximises the sum of the particle 
momenta projections in the two-dimensional plane orthogonal to the thrust axis: 



T M Q = max < V] | np h | > = V" | n M p h | = V] | p hy | . (2.3) 

nnT=0 

\ h ) h h 

The direction % we shall identify with the y-axis. Finally, for the thrust-minor we 
have 

T m Q = ^2\phx\ = K ont . (2.4) 

h 

We choose ny in such a way that the most energetic particle in the event has a 
positive ^-component and % in such a way that the second most energetic particle 
has a positive y-component. Hereafter we attribute ph to the right hemisphere Cr 
(left hemisphere Cl) if Phz > {phz < 0). Similarly ph is in the up hemisphere C\j 
(down hemisphere Co) if Phy > (ph y < 0). We shall also consider separately the 
right-hemisphere cumulative out-of-plane transverse momentum: 

It can be easily shown that from the definition (2.2-2.3) the kinematical constraints 
follow: 

^2pht = ^2pM = Q , ^2phx = ^2phx = v, (2.6) 

h<=R heL h£U heD 

where p t is the two-dimensional vector transversal to the z-axis. We introduce the 

— * — * 

three-dimensional vector Pi and the two-dimensional vector P 2t that define the "event 
plane" {z,y}: 

A = (P lx , P ly , P lz ) = (0, 0, TE) , P 2t = (P 2x , P 2y ) = (0, T M E) , Q = 2E. (2.7) 
From the definition of T and T M we have 

In what follows we study the distribution of events in the cumulative out-of-plane 
transverse momentum K out defined in (2.4). The integrated K out -distribution is 
defined as 



do[K Q 



dTdT M 



Q 5 J2 [ ^ m e(K out -f>^|U 3 ^ 

* h=i J \heR J \heu 



(2.9) 

where m denotes the number of final particles in an event. The last two delta- 
functions fix the event plane and the theta-function defines the observable. Analo- 
gously we define the right distribution by restricting the sum over particles in the 
theta-function of (2.9) to those belonging to the right hemisphere (see (2.5)). 
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Figure 1: The three Born configurations £3,^2 andCi for T = 0.75 and Tm = 0.48 ordered 
according to decreasing probability. The thrust T and thrust major Tm are along the z- 
and y-axis respectively. The four regions Ce in the phase space are indicated. 



2.1 Jet momenta and hard parton recoil 

At the parton level the events in the region (2.1) can be treated as three-jet events 
generated by a system of energetic quark, antiquark and a gluon accompanied by an 
ensemble of soft partons. 

At the Born level, O (a s ), in the absence of accompanying radiation the 3-parton 
system is truly planar, T m = 0. The kinematical configuration of q, q and g treated 
as massless partons is then uniquely fixed by the values of T and Tm- Denoting by 
Pi,p2 and P3 the energy ordered parton momenta, P±o > P20 > -P30! we have Pi 
lying along the thrust axis. The vectors Pi and P 2t = — P 3t are given by the event 
plane momenta defined in (2.7). 

There are various kinematical configurations of the Born system. The three par- 
ton momenta P q , P q , P g can belong to three configurations Cs, namely (see Figure 1) 

(Pi,P 2 ,P 3 ) = (P g ,P q ,P q ) => C 1 

= (P q ,P g ,Pg) => C 2 (2.10) 

= (P q ,Pg,P g ) C 3 . 

Notice that the index S labelling the configuration Cs coincides with the index of 
the gluon momentum. (Interchanging the quark and the antiquark does not affect 
the accompanying radiation and the distributions under study.) For three massless 
partons the values of T and Tm are restricted to the region 

2 (1 ~ T) v / 2T^T <T M < VT^T . (2.11) 

(For the kinematics of the Born system momenta see Appendix A.) In the following 
we will analyse the distribution inside this region — the 3-jet region — in which 
three skeleton parton momenta P a can be reconstructed from the T and T M values. 

Beyond the Born approximation, in the presence of secondary gluon radiation, 
T m is no longer vanishing. The a>components of bremsstrahlung gluon momenta are 
logarithmically distributed over a broad range, so that (T m ) ~ o: s Tm ~ ct s T. 
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In the region (2.1) the PT expansion develops logarithmically enhanced con- 
tributions which need to be resummed in all orders. Our aim is to perform this 
resummation with SL accuracy. This means that we shall keep both DL, a s log 2 T m , 
and SL, ct s logT m , contributions to the exponent ("radiator", see below) while ne- 
glecting non-logarithmic corrections 0(a s ). With account of the running coupling 
effect, the DL and SL contributions formally expand into series of terms <log n+1 T m 
and a"log ra T m , respectively. 

It is important to stress that exponentiation of the SL correction a s log T m makes 
sense only if it is supported by the calculation of the non-logarithmic "coefficient 
function" . Indeed, a "cross-talk" between an O (a s ) correction to the overall nor- 
malization of the resummed distribution and the leading DL term, symbolically, 

s 2 

(1 + c • a s ) xexp {a s log 2 +s ■ a s \og} = 1 +. . .+ c • a s x a s log 2 +-^(a s log) 2 +. . . , 

(2.12) 

gives rise to the correction of the same order as the SL term squared. The coefficient 
c is not known analytically. It can be determined by comparing the O (a 2 ) term of 
the log-resummed expression with the result of a numerical calculation based on the 
exact a 2 matrix element [10]. 

We denote by p a the momenta of the three hard partons, q, q, g, that in general 
no longer lie in the event plane defined by the vectors (2.7). In the region (2.1) they 
differ from P a by "soft recoil parts" q a , 

Pa = P a + q a , (2.13) 
whose re-components contribute to K out together with the soft parton momenta kf 

K out = \qix\ + \q 2x \ + 1 9ac| + ^2 ^ ' ( 2 ' 14 ) 

To SL accuracy, as in the case of the broadening distribution [3], it suffices to keep 
the recoil momenta q a only in the phase space, in particular their contribution to K out 
and to the constraints defining the event plane (see (2.9)). At the same time, q a can 
be neglected in the radiation matrix element, that is, in the emission distributions we 
can substitute the skeleton momenta P a for actual parton momenta p a . We remark 
that, as in the case of jet broadening [12], for the study of non-perturbative power- 
suppressed corrections the approximation p a — > P a in the radiation matrix element 
is not valid [13]. 

In order to resum the PT series for (2.9) in the region (2.1) we need to use 
the factorization property of the soft radiation matrix element and to factorize the 
multi-parton phase space. We discuss these points in succession after giving a brief 
description of the kinematics. For more details see Appendix A. 
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2.2 Kinematics of soft or collinear emission 

Introducing the recoil momenta q a defined in (2.13), the right, left and total K out are 
given by 

= I Six | ' K ont = l«2*l + \q 3 x\ , K out = K^ t + K^ nt . 

R L 

(2.15) 

The event plane momenta (2.7) are given by 

A = Pi + Y k ^ P 2t = p 2t + Pit$(Piy) + Y kit • (2.16) 
R U 

In the region in which the emitted partons are soft or collinear the components of q a 
that accompany non-vanishing components of P a can be neglected in the calculation. 
Therefore we consider only those four components that vanish in the Born approx- 
imation namely, the y-component of the leading parton momentum, qi y = pi y , and 
the x-components q ax = p ax , a = 1,2, 3. Then, from the kinematical relations (2.16) 
and (2.5) we have 

+ Y ki y = <7i* + Y klx = ' 

R R ^ 217 ^ 

<l2x + qtx + Y klx = q3x + q ix + Y klx = > qf x = q lx $(±q ly ) . 

U D 

Notice that this implies the following kinematical relations 

{Qtx + Y k ™) = + Y kix ) = ~i q2x + Y kix ) = i q3x + Y ki *}- (2.18) 

ieCi ieCA «ec*2 «eC3 

2.3 Matrix element factorization 

At the parton level the distribution is given by the sum of the partial cross sections 

do-n = -,M 2 n - d$ n , (2.19) 

71/. 

where M n is the matrix element for the emission of n secondary partons off the qq, 
g system, and d$ n the corresponding phase-space factor, 

<**« = (2W^p a +^^-g)nW-n-^-, [dk} = ^. (2.20) 

As it is explained in detail in Appendix B.2.2 the contribution of collinear non soft 
emission is process independent and can be taken into account a posteriori introduc- 
ing in the computed distribution the hard part of the splitting functions. 
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Therefore the accompanying partons ki are assumed to be soft and their distri- 
bution can be treated as independent. Assembling together the contributions from 
the three configurations Cs defined in (2.10), the distribution can be presented as 



8=1 i 

where M (C$) is the Born qqg matrix element and W$ is the distribution of the soft 
gluon radiation off the hard three-parton antenna in the momentum configuration 
Cs- Here we discuss the real emission contribution to M n ; the virtual corrections will 
be accommodated later. 

For the configuration 5 = 3, for example, the squared Born matrix element reads 

2 C F a s x\ + x\ _ 2P a Q 

M ° (Cs) = 2vr (l- Xl )(l- X2 ) ' X « = -Q2- ' (2 ' 22) 

where P 3 is the gluon momentum and Pi,P2 the quark-antiquark momenta (see 
(2.10)). For this configuration the single soft gluon radiation pattern is given, at the 
one-loop level, by 

W 3 (k) = C F w 12 + (^13 + w 23 - W12) = ^ ^13 + w 2 3 - j-pW 12 ^j . (2.23) 

Here w a b is the standard two-parton antenna of the a6-dipole, which, within the 
normalization convention prescribed by (2.20), is given by 

a s (P a Pb) a s 

w ^ = = ■ (2 - 24) 

Here k t)0i b is the invariant gluon transverse momentum with respect to the hyper-plane 
defined by the P a , P5 momenta. 

The first term CfW\ 2 in (2.23) is the "Abelian" contribution describing soft gluon 
emission off the qq pair. The second term proportional to iV c is its "non-Abelian" 
counterpart that describes radiation off the hard gluon P 3 . Similar expressions for 
two other kinematical configurations (5 = 1,2) is straightforward to write down by 
properly adjusting the parton indices, 

N e ) 1 '« < M5) 

W 2 (k) = — {w 23 + W12 - Jp w i3 

To reach SL accuracy it is necessary to treat multi-parton emission at the two-loop 
level. This involves allowing secondary gluon to split into two gluons or into a qq 
pair. In principle, perturbative analysis of a system consisting of three hard partons 
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and two softer partons (with comparable energies) can be found in the literature 
(see, e.g., [14]). However, to the best of our knowledge, the most important feature 
of the result has never been explicitly stressed, namely, that the colour structure and 
geometrical properties of emission of a soft two-parton system off a 3-jet ensemble is 
identical to those for single gluon radiation. 

It can be shown [15] that after subtracting the uncorrelated radiation of two 
soft gluons, W s (k 1 ) W 3 (k 2 )i the correlated two-parton production is given by the 
expression 

wP(k ± , k 2 ) = C F w { S(h, k 2 ) + ^ (w®(k u k 2 ) + w%\k u k 2 ) - w^(h, fc 2 )) , 

(2.26) 

(2) 

where w ah is the standard distribution known from the two-loop analysis of 2-jet 
event shapes (the first term on the right-hand side of (2.26), see [16]). It describes 
decay into qq- or gg-pair of a virtual parent gluon radiated, in our case, by one of 
the three two-parton dipoles ab. Given this remarkably simple dipole structure, the 
analysis of the two-loop effects in 3-jet events reduces to the known 2-jet case. 

As shown in Ref. [16], for a sufficiently inclusive observable (and the i^out- 
distribution under interest in particular) the two-loop refinement results in (and 
reduces to) substituting the proper argument for the running coupling describing the 
intensity of the parent gluon emission, 

Wabik) = -y^— =>• w ab {k) + [ dk 1 dk 2 5(k - h - k 2 ) w^(k u k 2 ) ~ "^'"^ . 

(2.27) 

Here k t , a b is the invariant transverse momentum defined in (2.24) and a s is taken in 
the physical scheme [17]. It is worthwhile to notice that the arguments of the cou- 
pling are different for the three dipoles that participate in gluon radiation according 
to (2.23). 

We remark that at the level of the leading power-suppressed non-perturbative 
correction the role of two-loop effects is more involved as they give rise to the Milan 
factor [6,16,18]. 

2.4 Phase space factorisation 

We discuss here factorization of the phase space in the soft region, which is needed 
for resummation. The phase space factor reads 

dT n = d§ n ■ 5 3 (p!+J2 k-P^j 5 2 (^2t+Pit${piy) + ^kit-P^J 



(2.28) 



1TL) 
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where the factor D n takes care of the kinematical relations, 

A^Tr)^ 4 ^^^*;^ . 

(2.29) 

The first delta-function stands for the energy-momentum conservation, while the 
last two (three- and two-dimensional) delta-functions set the event plane. 

Now we single out from dT n small recoil components of hard partons momenta, 
Qiy — Piy an d the three q ax = p ax , which have to satisfy four event-plane constraints 
given in (2.17). Then, introducing the unity, 

1 = / dqiy Y[ dq ax ■ S n , 

J a 

where 

(2.30) 

we neglect hard parton recoils q a i in all but these four components and approximate 
the kinematical factor D n as follows 

Dn = J dq ly Y[ dq ax S n • D n ~ D ■ J dq ly JJ dq ax ■ S n . (2.31) 

a "a 

Here D is a trivial phase space factor which corresponds to the Born three-parton 
kinematics, 

A) = (2tt)V (^p a - Q^j 5 3 (p! - A) S 2 (p 2t - P M ) • 

We then have 

n 3 3 ,q 

dT n ~ T T\[dki] dh n , dh n = dq ly TT dq ax ■ S n , r o = T f , n »n r - A) , 

- LJ : - LJ r {2,-KyZha 

i=l a=l a=l v ' 

(2.32) 

where To is the Born phase space given in Appendix A. 

Neglecting bremsstrahlung in the D -factor proves to be legitimate: it can be 
shown that to achieve SL accuracy it suffices to take care of accompanying parton 
momenta in the ^-factor (2.30) and in the observable itself. 

Finally, in order to "exponentiate" the multiple radiation we need to factorize 
dependence on individual secondary parton momenta contained in the delta-functions 
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in S n and in the i^ out -observable. This is achieved in a standard way by means of 
Mellin and Fourier representations. In the following we apply this procedure to the 
integrated i^ out -distribution defined in (2.9). 

In what follows we shall separately consider K ont accumulated in the right hemi- 
sphere and the total K out of the event. We start from a simpler case of the K out - 
distribution in the right (one-jet) hemisphere. The total i^ out -distributions will be 
considered in Sec. 4. 



3. Right i^ out -distribution 

Consider the differential three-jet cross section with given T and Tm and with accu- 
mulated out-of-event-plane momentum in the right hemisphere smaller than a given 
K ut- Using the soft-factorization formula (2.21) we can write the cross section (2.9) 
for small K ont in the form 



da R {K out 
dTdT M 



3 ^ (0) 



(3.1) 



where af^ is the three-jet differential Born cross section for the parton configuration 
C s defined in (2.10), 

r M 2 (c 5 ) , 



dTdTM 

calculated in Appendix A. The accompanying radiation factor E, a function of K out , 
T and Tm, reads 



H R (K ont ) = j dh n d ^out-l qi x I lj 



(3.2) 



We recall that within the adopted PT accuracy the hard parton momenta p a entering 
into the soft gluon distribution factor Ws(k) can be approximated by the event plane 
vectors P a so that integration over the recoil variables q a = p a — P a can be easily 
performed. Since the observable involves only momenta in the right hemisphere, 
the recoil momentum components in the left hemisphere q 2x and q 3x can be freely 
integrated out with use of the last two delta-functions in (2.30) and leave no trace 
in the distribution under consideration. Then, a non-trivial dependence on q ly (via 
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qf x )) also disappears, and the (^-integration trivializes as well. We are left with a 
single delta-function in (2.30): 

H R (K out ) = j dq lx § ^ftT ut-| Qix \~y~] 1 hx |j 6 ^qi x + k ix ^j . 

To factorize the dependence on secondary parton momenta we use the Mellin rep- 
resentation for the theta-function in K out and the Fourier representation for the 
delta- function in qi x : 

•& (k^x - 1 q lx | ~y] 1 k ix |j = J ^-exp jz/ • ^T out -| g lie |- | fc ix |j J (3.3) 

5 ^/ix+ ^ = " y ^~ ex P ' ^9ix+ ki *j | ' ( 3 - 4 ) 

where the contour C in (3.3) runs parallel to the imaginary axis at Reaiz/ > 0. Defining 
the Fourier integral (3.4) we have extracted, for the sake of convenience, the factor 
v as if it were a real parameter: for a complex value of v it implies rotating the 
/3-contour by — Argz/, | Argi/ | < 7r/2. (Then, using the analytic continuation, the 
/9-integral can be transformed to run along the real axis.) 
Integrating over q lx we obtain 

H R (K *) - f ^ C vK ™* 1^ — Tl e -v(\kix\+i/3k ix ) (or) 

The limit of small K out (and thus of small q\ x ) corresponds to the region of large 
values of the conjugate variables, i.e. v and z//3 respectively. Therefore in what 
follows we will concentrate on the limit of large v and neglect the contributions of 
the order of u^ 1 which correspond to O (K out ) corrections to the distribution. At the 
same time, by examining (3.5) it is easy to see that the characteristic values of the 
rescaled Fourier variable (3 are of the order of unity. 
Substituting (3.5) into (3.2) we get 

where we have introduced the "radiator" 

Kf(v,P) = [ [dk] W s (k) [1 - e -"(lfel+i^) ] . (3.7) 

J R 

Here the unity in the square brackets has been included to account for the virtual 
corrections. 
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For example, for the most probable jet configuration 5 = 3 (with gluon the least 
energetic parton), the soft distribution (2.23) gives 

n 3=Y ( r ?3 + ^ - ^r? 2 ) , r* = J[dk\ w ab (k) [ 1 - e -(l*.l+tf*.) ] . 

(3.8) 

With the DL accuracy only gluons collinear to Pi (the thrust axis direction) con- 
tribute to the right-hemisphere K out . Therefore the (identical) DL contributions are 
contained in rf| and rf|. According to (3.8), they combine into the expression pro- 
portional to Cp — the colour charge of the quark P 1 . Single logarithmic corrections 
to these two dipoles, as well as SL contribution of the third dipole, r||, are calculated 
in Appendix B.4. Here we report the result: 

r£ = r(p, Q 2 ) + F la [ Q , r* = p 23 f > (3 . 9) 

Jl/P k x 7T ~ J 1/p k x 7T 

where the variable fi originates from an approximate evaluation of the characteristic 
momentum integral which is explained in Appendix B.6, 

[ 1 - e^l+tf*") ] -> e(\h x \ - , /Z = P/x , /i = v / l+^, (3.10) 

with 



u = ue lE . (3.11) 

In (3.9) r(/Z, Q 2 ) is the DL function 



f Q 'dk x aA2K) Q' 2 

rM ) = A /P fc -- ^ IfeT' (3 ' 12) 



and the factors F a t = F a b(T, Tm) are independent of the integration variables z/ 
and /3. The origin of an essential subleading correction embodied into the precise 
argument of the running coupling, a s (2k x ), is explained in Appendix B.5. 

Combining these results we obtain the radiators, evaluated with SL accuracy, for 
each of the three kinematical jet configurations, 

Kf(j£) = C s r(/i, Q 2 ) . (3.13) 

Here Cg is the colour charge of the hard parton along the thrust axis, that is C$ = C F 
for 5 = 2, 3 and Ci = N c . SL corrections in (3.13) were absorbed into the definition 
of the hard scale. The corresponding scale Qs depends on the event configuration 
5 and, though the functions F ab (T,T M ), on the event kinematics. These scales are 
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given by 

N c XnQl = N c ln(Q 2 e-^) + ^ (F 12 + F 13 - ^F 23 ) , 

TV 1 
C F InQ 2 = C F ln(g 2 e-^) + (F 12 + F 23 - ^F 13 ) , (3.14) 

C F InQ 2 = Cf ln(Q 2 e-i) + ^ (F 13 + F 23 - -^F 12 ) . 

They also include the factors e -3 / 4 and e _/3()//4Arc coming from SL corrections due 
to the "hard" parts of the quark and gluon splitting functions. As we shall see 
below in Sec. 4, for the case of the total if out -distribution the corresponding scales 
are simply related with geometry of the 3-jet event. At the same time, the T- 
and T M -dependence of the scales entering into one-hemisphere distribution does not 
have a simple geometrical interpretation. This is due to the fact that the kinematical 
constraint restricting the observable to a single hemisphere is foreign to the structure 
of the soft-gluon radiation pattern. 

The radiator (3.13) depends on the Mellin-Fourier moments v and f3 only via the 
variable p. The /^-dependence can be further simplified at SL accuracy by expanding 
the radiator for large z/, 

nfip) = nf{v) + C s r>{V) In y^YW 2 , 

where 

r'(u) = ^\n%, k = l/P. (3.15) 

7T K x 

Since r' constitutes a SL correction, we can use Q as a common scale (neglecting 
O (a s ) mismatch) and omit the factor 2 in the running coupling argument (as pro- 
ducing a negligible correction O (a 2 logz/)) in such subleading terms. 

Thus, the right-hemisphere i^ out -distribution to SL accuracy takes the form 



Integration over the Mellin variable v can be performed by steepest descent or by 
the operator method introduced in [12]. We follow the last method which exploits 
the following identities: 



m = /(e-*).(p)- 

m+z)' 



z=0 



J 2^ e "° Ut = ' Koat = e ^ EKout 
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The following approximation is valid: 



p -R(e- d - ) <p( p -d x \ ( K outY 



z=0 

where we have neglected relative corrections of the order R"(x) = xd x R'(x) = O (a s ) 
and assumed that T is a smooth function of R! which is true for F^v) in (3.16), 
dlnF/dr] = 0(1). 

Applying these relations to our distribution in (3.16) we derive the final answer 
to SL accuracy, 

Sf (K ont ) = e^A) . v( fl^fL ^ , ^ut = e--X out . (3.17) 

T (1 + C s r f {K ont )) 

The first and second factors resum DL and SL contributions, respectively. We remark 
that the precise argument K out is essential to keep in the first factor, while substi- 
tuting K out by, say, K out in the second factor would amount to a negligible O (a s ) 
correction. Comparing the //-integrand in (3.16) with the final result we conclude that 
the SL factor T -1 in (3.17) accounts for a mismatch between the Mellin-conjugated 
v and -Kout values. It can be looked upon as a next-to-leading order prefactor of the 
WKB (steepest descent) approximation. 

4. Total if out -distribution 

The calculation is similar to the previous case except that now radiation in all four 
quadrants contributes to K ov &. The -fT out -integrated distribution (2.9) for small K out 
is given by 

da(K out ) _ x ^ daf ] T 
dTdT M ~2^dTdT M 5( ' outJ ' 

Cg 

r i n 

Zj(K out ) = J Y.~\ II W WW ' # T (^out) , (4.1) 

n i 

H T (K out ) = J dh n d ^-ftTout-y^ \q a x\-^2\ hx |j • 

To factorize the soft momenta in H T we proceed as before by using Mellin and Fourier 
representation for the theta- and delta-functions. Again we denote by v the variable 
conjugate to K out and study the region \v\Q 3> 1. We rescale by v the variables 
conjugate to the soft recoil variables q ax (a = 1, 2, 3) and q\ y in dh n to arrive at 



dh n JJe-"' 9 -"' JJe""!* 4 -! = 

a=l i=l 

/!£ (n^r) / (' 9 ^){n^(A: l )n M i3(A;,)n^^)n M 3^)|. 



(4.2) 
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Here we have introduced the probing functions u a (k) for each of the quadrants Cf 

u 12 (k) = «(/9i2,7) » wi 3 (fc) = m(/9i3, -7) , «2(fc) = 0) , « 3 (A;) = u(p 3 , 0) , 

(4.3) 

where p 12 = Pi + P2, P13 = Pi + P3 and the "source function" u is 

u(Pn) = exp fcia, I +i/3/c ia; + i7|A; i2/ |)}. (4.4) 

The function 7) in (4.2) is the result of integrations over the recoil momenta, 
namely, gi y (conjugate to 7) and three g aa; (conjugate to /5 n ): 

= T+M i + M + ^tT7 + TTph ii + d ' (4 ' 5) 

Notice that the large variable v enters only in the sources. 

Now we are in a position to resum multiple accompanying radiation. The result 
reads 

Es(A '°"') = / 2^ e " ou ' / % n f 7) ■ . <«> 

a=l 

with the radiator given by 

f 4 
P, 7) = / [dfc] [l-£ • ( 4 -7) 

As before, the unity in the square brackets has been included to account for the 
virtual correction contribution. Here Qg(k) is the support function for a parton k 
emitted in the quadrant Ci, and we have denoted U\ = u± 2 and 114 = u± 3 . 

In Appendix B the radiators 1Z$ are evaluated with SL accuracy and we obtain 

= r(j2 2 ,Ql) + Ci S) r (jl 3 ,Ql) + ^- J [r(jl 12 ,Ql) + r(j2 13 ,Ql)} , 

(4.8) 

where r is the DL function defined in (3.12). Some comments on the colour charges, 
the scales and the various /i-variables are in order. 

• The radiator consists of three "independent" contributions from the radiation 
off each of three hard partons. Here 

C^ = N C ; C® = C F , fora^fe, (4.9) 

is the colour charge of the parton P a . 
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• The hardness scale Q\ in each term has a simple structure: it is determined by 
the invariant transverse momentum of the hard parton a with respect to the 
dipole be: 

Q2 _ Pta -g a 2 _ 2 (PbPg)(PgPc) ^ 

This makes Q a depending on the event geometry. The scale Q\ also includes 
an additional, geometry-independent, factor e~ 9a depending on the nature of 
the parton a. This factor takes into account a SL correction due to hard parton 
splitting, with g a = 3/2 for a quark (a = 1,2 in (4.8)) and g a = (3 /2N c for a 
gluon. Due to this factor the scales depend on the configuration 5. 

Recall that the derivation of the exact form of the hard scales required a SL 
analysis which takes a due care of the inter-jet regions where the soft distribu- 
tion does not have collinear singularities. 

As we have seen in the previous section, the hard scales for the right dis- 
tribution (3.14) do not have such a simple interpretation because selecting 
one hemisphere is unnatural for the soft radiation pattern. As shown in Ap- 
pendices B.2.1 and B.2.2, the boundary effects due to radiation at 90° which 
complicate the scales, cancel with SL accuracy in the total distribution. 

• The various functions Ji originate from the following substitutions in the inte- 
gral over k x 

l-u a ^#(k x -l/ji a ), Jia = v V~l + PI = v ■ »a , a = 2, 3 (4.11) 
1 - u 12 -> -&[k x - I//Z12) , Ji\2 = v \J(l -i'jy) 2 + Pl 2 = y-^\2- (4.12) 
1 - u 13 -> d{k x - 1//Z 13 ) , fi 13 = v sj (I + i-iy) 2 + (3 2 u = v ■ // 13 , (4.13) 

with v given in (3.11). As in the case of the right radiator, such substitution 
is valid within SL accuracy. The 7-dependence enters only in the two Ji\ a 
parameters which are associated with the gluon emission off the parton P 1 in 
the right hemisphere. This is in agreement with the fact that the variable 7 
is conjugate to q\ y . Opposite signs of 7 in p, i2 and fti 3 reflect the fact that 
the recoil momentum q ly is positive (negative) in the right-up (right-down) 
quadrant. 

All v-, f3- and 7-dependence is contained in the parameters p,. 

As before, we can simplify the (3- and 7-dependence by expanding the various terms 
to SL accuracy for large v. We can write 

n s = Rs(is)+r'(v)-Ss((3, 1 ), (4.14) 



19 



where R$ is the DL contribution, given by the sum of three standard antenna terms, 
r' is the SL function defined in (3.15), and the coefficient Ss carries all the f3 a and 7 
dependence. Rs and Ss are given by 

Rs(P) = C? r(P, Q\) + C? } r(P, Q\) + cf r(P, Q*) , 

s 6 {p,i) = d!P ln^+cf 1^3+^ jH j^MMvWv)], (415) 

with the charges given in (4.9). In detail the DL terms for the three configurations 



are 



R 3 (P) = C F r(u, Q\) + C F r(P, Q|) + iV c r (P, Q*) ? ( 4 . 16 ) 
i? 2 (p) = C F r (P, Q?) + iV c r(P, Q') + Cp r (P, Q*) ? (4.17) 
(P) = iV c r(P, Qi) + C F r(P, <$) + C F r (P, Q\) . (4.18) 

There are two sources of SL corrections in (4.14). The first is due to different hard 
scales in three terms r(P, Qi) in (4.16)-(4.18), which depend on the geometry of the 
three-jet event, that is on the values of T and T M . The second is the contribution 
proportional to r' given by the sum of In /x-terms which depend on (5 a and 7. These 
contributions are specific for a three-jet topology. At the same time, within the SL 
accuracy this correction is insensitive to the details of the event geometry, that is to 
the values of T, T M ~ 1. 

In conclusion, the total if out -distribution, to SL accuracy, can be expressed by 
the following Mellin integral: 

S,(^ out ) = J ^~ e» K ^ e~ R ^ ■ F s {y) , 
where the SL pref actor J 7 sip) is given by 

= P- fl r — • W> 1) e- r ' (P)S ^. (4.19) 
J —00 ^ 1 J —00 71 



a=l 



For example, the explicit expression for JF 3 reads 



a=l 

x exp |-C F r'(P) jT In [^(j)^)]}} • 

The function JF 5 is analysed in Appendix C. 

In the first order in a s \nu the factors T become (see Appendix C.l) 

Ts = 1 - (3C F + JV C ) In 2 • r'(p) , 5 = 2, 3 , 
^1 = 1- (2C F + 2N C ) In 2 • r'(p) . 
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Different weights of the quark and gluon colour factors for the two cases have a simple 
explanation. Due to the kinematics of parton recoil (see (2.17)) contribution of the 
gluon radiation off the most energetic (right-hemisphere) parton Pi is twice that off 
the left-hemisphere partons P 2 and P 3 . As a result, the SL correction is proportional 
to 2Cf + Cf + N c when Pi is a quark/antiquark (5 = 2, 3) and to 2N C + Cf + Cf 
when it is a gluon (5 = 1). 

Integration over the Mellin variable v can be done as before and we obtain 

S 5 (^ out ) c e-*(*£) ■ r( f| ( ^l t)) , K ont = e-^K out , (4.20) 

where R' is the logarithmic derivative of Rs- To SL accuracy we can take 

R'(K^ t ) = (2C F + N c ) ■ r'(K^ t ) , (4.21) 

which is the same function for all configurations S. This is possible since, r' being a 
SL function, the difference between the hard scales can be neglected at the level of 
the next-to-next-to-leading O (a s ) correction. 

4.1 Radiators in the quasi-2-jet limit 

In the 3-jet kinematics we have been considering, Tu ^ T < 1, jet energies are 
comparable and relative angles between jets are large. In these circumstances three 
basic scales in (4.10) are of the same order, Qi ~ Q 2 ~ Qs < Q. Still, keeping precise 
scales in (4.16)-(4.18) is essential, since their deviation from the overall hardness 
parameter Q in the DL radiator function r(P, Q%) produced a SL correction S^R oc 
r'(P, Q) ■ la(Q a /Q) = O (a s logz/). At the same time, since the next order expansion 
terms are negligible, 8^R oc r"(P, Q) ■ In 2 (Q a /Q) = O (a s ), it is perfectly legitimate 
to present the answer in a form different from (4.16)-(4.18). 

For example, expanding the scales in the first two (quark) terms in R 3 around 
q 2 = 2PiP 2 we obtain, instead of (4.16) 

- 2C F r (^e-§) + N c r(Ql) , Qi = 2 ^^\ -^ , ( 4 .22) 

where we have suppressed the first argument v and used In 2 pfp 2 = ~~ m 2pfp 2 ^° 
cancel the linear expansion terms, see (4.10). Mismatch between the radiator in 
(4.16) and the right-hand side of (4.22) is 

J R 3 (4.22)- J R 3 (4.16) = O (a s \n 2 . 

Being equivalent within the adopted accuracy, the representations (4.16) and (4.22) 
start to significantly differ, however, when the jet configuration becomes 2-jet-like. 
Indeed, when the gluon jet P 3 becomes relatively soft and/or collinear to the quark 
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P 2 , we should expect the answer to correspond to the gg-dominated radiation pattern, 
and the gluon contribution to disappear. 

The latter representation (4.22) correctly describes this situation: the quark- 
antiquark antenna contribution with the hardness scale equal to the invariant squared 
mass of the qq pair, q 2 = 2PiP 2 , takes on the job, while radiation off the gluon 
vanishes with decrease of the gluon transverse momentum, Ql ~ p 2 3 <C q 2 ~ Q 2 . At 
the same time, the first representation goes hay-wire. For example, in the collinear 
limit, P3P2 <C P3P1 ^ P\Pii one of the quark scales, Q\, vanishes (together with 
Q 2 ) while the other formally goes to infinity, Q\^> Q 2 . The geometric mean of the 
quark scales, 2PiP 2 = \Jp\ t v\f employed in (4.22) cures this unphysical behaviour. 

For two other jet configurations, with the gluon having an intermediate or the 
largest of the three parton energies, analogous representations of the radiators read 

R 2 ^2C F r ( + N c r{Q 2 ) , Q\ = 2 ( P f*Yf*) e -& ; (4.23) 



Ri ^N c 



(P1P3) 

I -±-± e -U+4N- c )\ + r I _L± e -(4 + 4iv^)J _ _ r (.4.24) 



These expressions, contrary to the original ones (4.17) and (4.18), survive the limit of 
the left-hemisphere jets becoming quasi-collinear, P2P3 <C P1P3 < -Pi Pi- Indeed, the 
configuration 5 = 2 is then identical, colour-wise, to 5 = 3, and (4.23) is dominated 
by the qq two-jet contribution, while the non-Abelian part vanishes with the gluon 
transverse momentum, Q\ ~ p\ t <C Q 2 . A rare but interesting configuration 5 — 1, 
where the gluon Pi in the right hemisphere is balanced by a quasi-collinear qq pair 
in the left hemisphere, see (4.23), corresponds to the gluon-gluon system: radiation 
off the colour- octet qq pair is dominated by large- angle coherent bremsstrahlung 
proportional to the gluon charge, N c , with the colour-suppressed 1/N C correction 
term vanishing in the collinear limit. 

Thus, the radiators in the form of (4.16)-(4.18) are inapplicable in the quasi- 
two-jet kinematics. An attentive reader could have noticed that the modified expres- 
sions (4.22)-(4.23), though better behaved, cannot pretend to uniformly preserve 
the desired SL accuracy. Accommodating correctly SL effects due to large-angle soft 
bremsstrahlung, these expressions fail, however, to properly account for "hard" SL 
corrections in the collinear limit. 

For example, the right-hand side of (4.22) for the most natural configuration 
5 = 3 has a perfect sq/t-gluon limit, when 2(PiP 2 ) = Q 2 becomes the proper 2- 
jet scale. At the same time, when P 3 remains energetic but collinear, E 3 < E 2 , 
(P 2 P 3 ) — > 0, the quark-jet scale in (4.22) remains smaller than the total annihilation 
energy, 2(PiP 2 ) < 2 (Pi ■ (P 2 + P 3 )) = Q 2 . The mismatch amounts to a SL correction 
0(r'-\nx 2 ). 
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5. Discussion and conclusions 

In this paper we performed the all-order perturbative analysis of the out-of-plane 
transverse momentum distributions in three-jet e + e~ annihilation events. We con- 
sidered the total i^out of the event and K^ ut accumulated in one event-hemisphere 
which contains the most energetic of three jets. The perturbative expression for the 
integrated -K" ut-distribution in the right (one-jet) hemisphere is given in (3.12)-(3.17) 
and (B.35). The i^ out -total distribution is determined by expressions (4.8)-(4.21). 
These answers resum all double- and single-logarithmic contributions to the exponent 
(the so-called "radiator" ) of the distribution in the standard Mellin- Fourier param- 
eter space. DL and SL contributions to the radiator can be formally represented, 
respectively, as series a™ \og n+1 (Q / K out ) and a™ log n (Q/-Kout) (n > 1). 

Throughout the analysis, we systematically neglected next-to-next-to-leading 
terms in the radiator (radiator series a™ log n - 1 ln(Q/K out ), n > 1) as well as non- 
exponentiating corrections of the relative order 0(a B ). The latter belong to the 
coefficient function 1 + ca s + O (ctg). Its first coefficient, c, is analysed in [10] where 
the order expansion of the approximate resummed cross is compared with nu- 
merical calculation based on the exact matrix element. Matching the resummed 
logarithmic expressions with the exact result is necessary for justifying exponentia- 
tion of the next-to-leading SL terms a s \n(Q/K out ). 

Accompanying gluon radiation pattern follows the colour topology of the under- 
lying parton antenna. The corresponding patterns were known, to SL accuracy, for 
two-parton — qq and gg — sources (quark and gluon form factors, 2-jet shapes). In 
the present paper we presented the first such analysis for 3-parton ensembles. 

The ggg-initiated events possess a rich colour structure which determines sec- 
ondary parton flows and makes them event-geometry-dependent. We have demon- 
strated that after taking into account SL effects due to inter-jet gluon flows, the 
result can be cast as a sum of quark- ant iquark (Cp) and gluon (N c ) contributions 
with the proper hardness scales depending on T and T M . 

A significant dependence on event geometry is the key feature which singles out 
the K out distribution among other e + e~ event shape observables. 

Soft gluon field components with small transverse momenta ( "gluers" with k t > 
^■qcd) are believed to be responsible for hadronization. This belief, known as hypoth- 
esis of the local parton-hadron duality, has been verified in a number of experimental 
studies of various features of multiple hadroproduction in hard processes (for review 
see [19]). Moreover, in recent years a new theoretical techniques have been devel- 
oped for triggering and quantifying genuine confinement effects by studying gluer 
radiation. 

A rich dependence of gluon (and, therefore, gluer) radiation on the colour topol- 
ogy makes 3-jet observables, and K out in particular, an interesting field for the study 
of non-perturbative effects. A separate publication [13] will be devoted to analysis of 
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the leading non-perturbative power-suppressed corrections to 3-jet observables. The 
K OVLt distribution, similar to the case of 2-jet Broadening observables [12], will be 
shown to exhibit log K ont -enhanced 1/Q contribution with the magnitude depend- 
ing on 3-jet geometry, that is on T and T M values. 

To access this interesting physics experimental studies of K out should be carried 
out specifically for events with moderate values of 1 — T, away from the 2-jet region. 

Two sorts of experimental studies of these predictions can be envisaged. The 
most straightforward comparison calls for experimental identification of the gluon jet. 
Gluon tagging, however, is unnecessary for the study of the total i^ out -distribution 
for genuine 3-jet events: the prediction given by the sum of three jet configurations 
corresponding to given T and T M , weighted with the proper 3-jet Born cross section 
factors bears practically as much information. Tagging brings in more information 
when a single-jet (right-jet) /Tout-distribution is studied. In this case essentially 
different i^ out -spectra will be seen depending on whether the right-hemisphere parton 
is a quark or a gluon. 
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A. Kinematics 



A.l Kinematics of qqg jets 

We consider the quark, antiquark and gluon with skeleton momenta P a in the centre- 
of-mass frame, J2l=i P a = Q = (Q, 0, 0, 0), 

Pi = %,0,0,ti), P 2 = E(x 2 ,0,T M ,-t 2 ), P 3 = E(x 3 ,0,-T M ,-t 3 ). (A.l) 

We have 2E = Q, x a = 2(P a Q) / Q 2 , x 1 + x 2 + x 3 = 2 and ti = t 2 + t 3 . Assuming 
x 2 > x 3 , we obtain 
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—P, (A.2) 
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where 



P=\ll-^ f <l. (A.3) 
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Thrust-major is confined to the kinematical region 

2(1 ~ T) v / 2T^T <T M < VT^, (A.4) 

where the upper limit comes from reality of p, and the lower limit comes from re- 
quiring x\ > x 2 . It is straightforward to show that in this kinematical region t 3 > 0, 
i.e. that both P 2 and P 3 lie in the left hemisphere. 

We introduce Sudakov variables based on two light-light vectors aligned with the 
thrust axis, 

P = E(l, 0,0,1) , P = E(1, 0,0,-1). (A.5) 

The jet momenta have the following Sudakov decomposition: 

P!=TP, P 2 = A 2 P + B 2 P + P t , P 3 = A 3 P + B 3 P-P t , (A.6) 

where the transverse momentum with respect to the thrust axis is P t = (0, 0, ET M , 0). 
In terms of T and T M the longitudinal Sudakov momentum components are 

A 2 = ^(l-p), B 2 = \{l + p), 

l-T I (A ' 7) 

As = —^-(l + p), B 3 = -(l-p). 

Since P a with a = 2, 3 belong to the left hemisphere, we have A a < B a . 

Denoting by a the angle of P a with respect to the thrust axis we introduce the 
angular variables r a 

H = t an —- -A. - 1 +p 

Q 2 QA 2 T M f \ a\ 

r 3 Qs -Pt I-P (A ' 8) 

— = tan — = — — = — — — . 
Q 2 QA 3 T M 

We remark that r 3 is negative since P 3 lies in the third quadrant ("down" hemi- 
sphere). Since P 2 and P 3 are in the left hemisphere and | P 2z \ > | P 3z |, we have 
0<7r — 2 <6 3 — 7r< \n. Hence, 

^>^>1. 
Q Q 

The following relations hold: 

T 2 I r 3 I = , t 2 + I r 3 | = — . (A.9) 

1 — 1 ±M 

The minimal value of thrust-major, T M = 2(1 - T)/(2 - T) for a given T in (A.4), 
corresponds to the configuration with the softest parton momentum orthogonal to 
the thrust axis, 

_ J_ T2 J_ Jj73_[ ^ 

P ^ 2-T ' Q ^ l-T' Q ~^ 
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The maximal value T M = y/1 — T is achieved with a symmetric 2 3 pair: 

n t 2 \t 3 \ 1 



Q Q VT^T' 
We also introduce two-dimensional vectors in the transverse plane {x, y}, 

r 2 = £ =(0,75), T3 = ^ = (0,T 3 ). (A.10) 

Since Pi is along the thrust axis we have t\ = (0, 0). 
A. 2 Born cross sections 

The squared Born matrix element M (Cs) has a well-known expression in terms of 
the variables x a . For 5 = 3 (with P 3 the gluon momentum) we have 

MHO = g«.W) (T^fe) • (A ' U) 

The phase space factor T for the Born 3-parton system is 

The dipole invariant masses in terms of T, Tm and variables r, are 



T 
2" 

-.2 _ 



Q 2 2 = 2P X P 2 = -(1 + P )Q 2 = \TT M r 2 Q , 



Qi 3 = 2P X P 3 = -(1 - p)Q 2 = ITT m \t 3 \Q , (A.13) 
Q 2 23 = 2P 2 P 3 = (1 - T)Q 2 = -^T • 

The scales in (4.10), 

2 _ 2(P fe P a )(P a P c ) 
Pta " (P 6 P C ) ' 

- the invariant transverse momentum of the hard parton P a with respect to the 
frc-dipole — read 

2 Q 2 ( TT M \ 2 2 2 /l + p\ 2 2 n2n rr^f^-P^ 1 

Pti = ir\ — ™ , Pt2 = <2 (1 - T)[ — — , p t3 = Q (1 - T) 



2 V 1 - T J ' 7 V T M ) ^ ' \T M 

(A.14) 
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A. 3 Soft partons 

For secondary massless parton of momentum k we write the Sudakov representation 

k = aP + f3P + k t , «/5 = t4, (A.15) 

Q 

with k t the transverse momentum with respect to the thrust axis. The right- 
hemisphere condition imposes the restriction upon longitudinal parton components, 

h 

a > p =>- a > — . 

Probability of soft gluon emission off the P a , Pb dipole is described by the squared 
matrix element 

p - n - — ° 2( y; t)2 - 2 . (A.i6) 



B. Radiators 



Here we compute the radiators for the total K out distribution, TZs(u, (3, 7), to SL 
accuracy. The radiator is given by the combination of dipole contributions, 

r ab = I [dk] w ab (k) [1 - H , (B.l) 

with Qe(k) the support function restricting gluon momentum k to the quadrant Ce 
and the source functions ue defined in (4.3), (4.4). 

Expressing the parton phase space [dk] in terms of Sudakov variables and invok- 
ing the soft distribution w a b in (A. 16) we have 

= f Q dk f Kym dky f amda l as ^ a2 ^ a ~ ^ 2 } vn lea) 



(B.2) 



We recall that f a are jet transverse direction vectors, T\ — and r 2j 3 given in (A. 10). 

We have three integrations to perform. The result of the first two integrals, in 
a and k y , has the structure 

Jkyjda oc ^(ln^ + 0(kS), ^ = 0{1), (B.3) 

where we have omitted O (k x ) terms as producing non- logarithmic O (a s ) correc- 
tions upon /^-integration. (For the same reason we have set the lower limit of 
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a-integration, a > (k t /Q) 2 , to zero, since the integral converges.) Similar correction 
originates from the precise limit of the k x integration, \k x \ ^,Q. Therefore we were 
free to set it to Q in (B.2). 2 

At the same time, the exact limits of a and k y integrals do matter when they 
affect the scale of the logarithm, Q', in (B.3) which are responsible for the SL correc- 
tion to the radiator. The logarithmic terms in (B.3) originate from three "collinear 
regions" when the gluon transverse direction vector f = k t /a is close to one of f a . In 
the direction Pi it is the a-integration that produces a collinear logarithm. Noticing 
that P\ z = TQ/2, we then set a m to the maximal value of a kinematically allowed 
in the right-hemisphere jet, 

a m = T. (B.4) 

In the regions collinear to the left-hemisphere jets a = 2, 3, a and k y are linked by 
the condition f ~ f a . In this case the kinematical limit can be expressed in terms of 
the ^/-component: P<i y = \ Pz y \ = TmQ/2 translates into 

T , T M Q 

K ym = (B.5) 

In (B.2) we employed soft radiation probabilities w ab . However, in the collinear 
regions hard parton splitting should be accounted for which produces another im- 
portant SL corrections similar to those coming from the kinematical limits (B.4) and 
(B.5). These corrections will be taken care of in the end of Appendix B.2. 2. 

We now proceed with successive integrations. Since the sources u do not depends 
on a we first compute the a-integral. Then we compute the k y - and /^-integrals. 

B.l Integrating over a 

We calculate separately contributions from the right (R) and left (L) hemispheres, 
R ^<a<a m = T, L a < ^ , (B.6) 

and define 

i* = r „ ada {fa i 7b? , i\ = [ kt/Q s ada ^z n)2 . (b.7) 

J kt/Q (h- ar a ) 2 (k t - ar b ) 2 ' a Jo (k t - af a ) 2 (k t - af b ) 2 
The a-integration yields 
ada (f a - T b ) 2 



o (k — af a ) 2 {k — afb) 2 

ln^-W + f* + $L ( ar ctan + arctan ^zf] \ 

+ h 2 ab kl\2 (ar b -k y ) 2 + k 2 x \k x \\ \k x \ \k x \ J) 



(B.8) 



2 We remind the reader that all these O (a s ) uncertainties, as well as those we shall encounter 
below, are related to the coefficient function [10]. 
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Here h ab is a function of r a and r b which has the following geometrical meaning: 



h a b — 



r a + r b t an \®a + tan \Q h sin^(9 a + 9, 



T a -r b tan \6 a - tan \Q h sin±(0 a -0 6 ) 



(B.9) 



— * 

with 6 a the angle between the momentum P a and the thrust axis. We have h al = 1 
(a = 2, 3) and h 2 3 = p. 

In the limit of small aplanarity, | k x \ <C Q, for the dipoles involving the right- 
hemisphere jet Pi we derive 



a 2 T 2 

a 



I L - — 

la ~ V 



2 K< a + 

2 k? 



kl 
kl 



+ 



k,, 
k, 



( 



k x | \ 



^ (1 - 2-&(-T a )) - arctan 



•ctan , ° , + arctan — ^ 



For the left-hemisphere dipole we get 



1 



kl + p^kl 

i 23 — ;_ 2 , „2J,2 



P hA + £ln 



P ln ^3 + g 

2 k 2 + /c 2 . 



arctan 



«2 



+ arctan 



, (B.10) 
(B.H) 



«3 

k T . 



kl + p>k\ 



P i K 2 + kl 
2 «2 + A; 2 . 



k x I \ 



/ k 2 «3 & v 

I arctan — — r + arctan — — r + 2 arctan — — 



(B.12) 



Here 



K 2 > , K 3 < . 



B.2 Integrating over k y 

To obtain the /^-integrand we need to perform the /^-integration 

rQ a s (2\k x \) 



Tab 



dk n 



IT 



B ab (k x ) . 



(B.13) 



The origin of the running coupling argument in (B.13) is explained below in Ap- 
pendix B.5. 

To compute the functions B ab (k x ) we consider two terms B^ and B^ coming 
from the upper (U) and lower (D) hemispheres, each containing left- and right- 
hemisphere contributions, 

" A "'"' W '" RVfl-^] [1-1*12] }, (B.14) 



(10 



7T 



and 



5 



(£>) 



r° dk 

= / ^ Wb-ll-us] [l-uis]}- (B.15) 
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From (4.3) we have 

«i2 = u(/3i2,7)> "is = u(Pis, -7) , m 2 = m(/9 2 ,0), u 3 = u(/? 3 ,0), (B.16) 
where /3i a = f3\ + f3 a and 

«(/3, 7 ) = e -"[l fc -l+^ fc -+*Tl fc »l] . (B.17) 
Only the right-hemisphere sources u\ a depend on k y . 
B.2.1 Dipole 23 

Consider first the 23-dipole contributions B^ D \ The DL piece is obviously con- 
tained in the left-hemisphere piece I^- At the same time, the contributions involv- 
ing are subleading: their fc y -integrals are not enhanced by log | k x |. Therefore, 
with SL accuracy the accompanying [1 — w la ] factors in (B.14) and (B.15) can be 
simplified, set equal to [1 — u a ] and factored out: 

lg-[l-u a ]+lg-[l-u la ] - (4 + /i)-[l-«J. 

Then, the terms in and J^ 3 due to the boundary between the R- and L-hemispheres 
cancel in the sum, and, using that the sources are now /^-independent, we obtain 

B$P = [1-U2].C, = [l-u 3 ]-C; 

dk f Pl r| 2k k \ ( B - 18 ) 

In — 7T + — — r arctan ■ 



7T(P+ P 2A;2) ^ 2 \k x \ \ k x 

where k = k y (k = —k y ) in the up (down) hemisphere. We extract the logarithmic 
contribution coming from the arctan term in the region | k x | <C k <C K m = \TmQ 
(see (B.5)) and using 



Jo 



/ 7T(F+p 2 ^) k 

we finally arrive at 



2h ' dk arctan l^i = In ±±£ , (B.19) 



where we have used the relations 

T 2 M = (l-T)(l-p 2 ) and 

following from (A. 3) and (A. 8), respectively. 



■2 n ,,,, :2, and 2L = l±P 
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B.2.2 Dipoles 12 and 13 

The contributions la can be simplified in a similar way. The left-hemisphere contri- 
bution (B.ll) is accompanied by ^-independent source u a . The same source can be 
attributed, however, to the second term in (B.10) as well, since the /^-integration 
is non-logarithmic and produces a subleading contribution. After this simplification 
the i?/L-boundary terms cancel in the sum of I^ a and l[ a , and we arrive at 



dk„ ( 1 , aLr 2 , r _ , k„ 



7T , . k, 



k x 



-sgn(r a )+arctan 
— sgn(r a ) +arctan 



y 



k x 



[l-« 2 ](B.2l) 



As before, it is important to keep the dependence of the sources on the finite pa- 
rameters j3 and 7 only in the contributions enhanced by ln(Q/ \ k x |). Such factor 
(collinear logarithm due to the jet #1) is explicitly present in the first terms in 
(B.21), (B.22). Adding these contributions together and integrating over k y results 
in substituting the R-hemisphere sources by the average source, 

[i-«i] ln Q 2 



b I £-2 



where 



l-«i = -/ -1 dy 5 [(1 - U12) + (1 - u 13 )], y = I k y /k x I . 
Jo 1 + y 

Collinear enhancement factor due to the jet #2 originates from logarithmic k y - 
integration in the second term of (B.21) (second quadrant), due to the jet #3 — in 
the second term of (B.22) (third quadrant). Remaining finite pieces can be evaluated 
using (B.19). 

Setting the limits a m = T and K ym = T M Q/2 according to (B.4) and (B.5) and 
using (A. 8), we finally obtain 

Bu = < +D > = ^ In % + ^ % + ^TT~p I" ^ . (B.28) 

Z I rv X I Kj, Z I K x I K x I /Cx I oC^ 

where wo i n the last subleading term is a source whose /3-, 7-dependence can be 
chosen arbitrary. Using this freedom we can absorb the last term in (B.23) into 
rescaling of the first two namely, 

We observe that the hard scales in (B.20) and (B.24) have a simple geometrical 
interpretation. Indeed, 

(1 - T)Q 2 = 2P 2 P 3 = Qj 3 , TTM ' Ta|Q = 2P lPa = Ql , (B.25) 
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with Qat, the invariant dipole masses, see (A. 13). 

Our analysis was based up to now on the soft radiation matrix element, (A. 16). 
To fully take into account SL effects from the region of large secondary parton mo- 
menta, we have, in addition to fixing the upper limits a m and k ym , to consider also 
hard collinear parton splitting. Due to collinear factorization, these corrections are 
process-independent and can be easily taken into account by proper rescaling of jet 
hardness parameters. They amount to supplying the invariant dipole masses by the 
factors 

f | for a quark/antiquark , 
Q 2 ab -> Ql b .exp{-U9a + 9 b )}, 9a ={ I r , (B.26) 



for a gluon . 



We finally obtain 



[(1- M2 ) + (1- M3 )] Ql 3 e-2^^ 

23 = ^PM ln M 

Bia = + (!-««)] ^Qle-fo*** 



(B.27) 



2 I k x | Akl 
B.3 Radiator by assembling bits and pieces 

Now that the a- and ^-integrations have been performed, we are in a position to 
assemble the full radiators for three jet configurations: 

K3 = £ ( Bl3 + B . 23 _ -L Bl2 ) , (B.28) 

K, = £ ik fJmi^ ( Bu + fi23 - ^sl) , (B.29) 

Rl = £ ^mL± ( Bl2 + Bl3 _ j_ B ^ . (B .30) 

The answers assume a simple when expressed in terms of the invariant transverse 
momentum of the hard parton P a with respect to the frc-dipole 

2 _ QbaQac 
Pta — n 2 
^bc 

It is straightforward to verify that, to SL accuracy, the answer can be represented in 
a symmetric form as 



n s = f * (cf m^[i - aj + a? [i - «j + cf > m%i - „,]) , 

7-0 «x 2tt V hi hi hi J 

(B.31) 
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where the hard scales are given by: 

Q\ = ?fe-^ , Q\ = P fe~^ , Q\ = ^ e "» . (B.32) 
The rule for the colour factors in (B.31) is simple: 



C^ = N C - Cf = C F , fora^b, (B.33) 

and the hard-splitting rescaling factors e~ 9a are defined in (B.26). 

The universal representations (B.28)-(B.30) do not have a smooth 2-jet limit. 
In subsection 4.1 alternative formulae are presented which are equivalent to the 
previous ones in the 3-jet kinematics, Tm < T = (9 (1), but are better behaved when 
the system assumes a quasi-two-jet kinematics, T M <C T. 

B.4 Radiator for the right distribution 

In this case the source does not depend on k y , and we have 

J - Q 71 (B.34) 

B? b (k x ) = [1 - e-KW-M*.)] ■ ^ Ml J* , 

J-oc ^ 

where the functions are given in (B.10) and (B.12). The DL contribution is 
contained in the first term of I^ a in (B.10). Extracting the large logarithm ln(Q 2 /k x ) 
which is embodied into the DL function r(/2, Q 2 ) in (3.9), we calculate the geometry- 
dependent SL correction factors denoted there by F a t,. 

The k y integrals are convergent so we have set the upper limit K ym — > oo. 
Introducing the ratio of momenta t = k y / \ k x | we obtain the following expressions: 

+ t (| sgn(r a ) - arctan«)) } , ^ 

— t (arctanfftf,) + arctan^g))} . 
where k' = = VT+ t 2 ^ — t. We remind the reader that convergence of the 

a \k x \ Q ° 

integrals is assured by r 2 > 1, r 3 < — 1. 
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B.5 Running coupling 



The two-loop analysis [15] entails that the argument of the running coupling for 
the dipole distribution r a b in (B.l) is given by k tja b, the invariant gluon transverse 
momentum with respect to the afr-dipole (see(2.24)). We show here that, to SL 
accuracy, we can the scale at the value 2\k x \ instead, as has been stated in (B.13). 
This effective value is obtained after integrating a s (k tab ) over a and k y . To show 
this we first observe that we need to control the precise argument of a s only in the 
DL contributions which originate from the phase space regions where the gluon is 
collinear to one of the three hard partons P a . 

Consider first the case of the contribution of r\ a from the right hemisphere. 
Here the soft gluon is close to the thrust axis, P±, so that the invariant transverse 
momentum reduces to the usual 2-dimensional momentum, k 2 lo ~ k 2 = k 2 + ky, 
which is a-independent. This allows us to perform the a-integration and obtain I^ a 
in (B.10). 

Now, to determine the effective scale of a s it suffices to consider fc^-integral of 
the leading piece of I^ a proportional to ha.(Q 2 /k^) and integrate over k y . We have an 
integral of the type 

.^^'kL^TT) ^^- (R36) 

with k x positive, y = k y /k x . For small k x the upper limit of the y-integral is large, 
y m = K ym /k x = T M Q/2k x ^> 1, and can set infinite since the integral converges. 
Expanding the coupling to the first order, we get 



K-A = jT -^ iy + — (a s (k x )-^a 2 s (k x )\n(l + y 2 )^ + O (« s 3 ) 



a s (k x ) - ^-a 2 s (k x )\nA\ - r\ s (2/,-.,) . 



(B.37) 



The regions collinear to P 2 or P 3 seems more complicated since here k ty ab depend 
both on a and k y . However, a similar analysis can be carried out in terms of the 
Sudakov variables with P, P aligned with the emitting parton, we obtain the same 
result for the argument of a s . 

B.6 The sources and the k x integrals 

Here we prove the substitution rule (4.11). The general structure of the radiators is 
D= ^\n^f[l-u], u = e- vk * cos{vpk x )e lv ™ k * , (B.38) 

JO kx 
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with Q' — O (Q) a hard scale. This we write as 



D _ r dk x ^q' | A A _ -^H I /p dk x (k x v r Q dk x n x y 

Ji/p, kx k x de \J k x \Q' J Jo k x \Q' J 



u 



e=0 

(B.39) 



and optimise the choice of p such that A = O (1), i.e. it does not contain a lnz/- 
enhancement. Since vQ ^> 1, the second integral containing the exponential source 
function, u oc exp(— vk x \ can be safely extended to infinity. 

Evaluating the expression in the curly brackets up to O (e), we obtain 

{ } = e- 1 [(jlQT - r(l + e) (vQT -Hi 1 - i-W ~ iP)" + i 1 ~ <7V + tf) _e )" 
= - (vQT (in £ - IE ~ \ In [(1 - iiy) 2 + P 2 ] + O (e)) . 

(B.40) 

Taking the e-derivative in (B.39) we obtain a large parameter \n(vQ') accompanied 
by the factor which we set equal to zero to optimize the choice of p: 

p = ve^y/il -iiy) 2 + I3 2 . (B.41) 

This means that, within SL accuracy, the source factor [1 —u] can be substituted by 

[1 _ u ] _ &( kx _ . (B.42) 

Performing this substitution in (B.31) we get the result reported in the text, see 
(4.14). 



C. Evaluation of T 

The expression for J-'g is rather complicated. Invoking (4.5) we split the 7-integral 
into two pieces, 

T = Fr + Fi, (C.l) 
namely the principal value and the 5(j) contributions, 

! - ± in5(-f) . (C.2) 



7 =F xe 7 

We have 



T _ r d A r d A n + T \ ,„ ^ 

J ~°° 7r(l + /?|) 1+ 2 C 2 ' J -co ^(l _|_ f]^y+2 C 3 r 7T 
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where the two integrands are 



c\ 5 K> 



1 + /% 1 + 



°^ (7( i +^ 2 +/W( i+7 ) 2+ ' 9 



2 ^3 



sin(Cf } r'Ai 



7 



Here 



1 Z" 00 
2Wo 1 



In 



(l+ 7 x) 2 + /3 

(l+7) 2 + /3; 



27T J 1 — X 

Where as before the colour factors are given by (B.33). 

C.l T in the first order 

In the O (r 1 ) approximation we have, for T r 

d(3 2 r d(3 3 



12 _ ln (l + ix) 2 +f3] 



(l+7) 2 + /5; : 



(C.4) 



(C.5) 



/oo 
- 
■oo 7T 

Using 



/■°° ^ 3 /- 00 i \ i 



l ■ 



Wo 7 



8~l + /3? ' 
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we find 
Using 



r'c?_ r°° d(3 2 r°° df3 3 r°° d^_ ( i i \ 

~J-ocAl + Pl) J oo W+W) J-oo 2^ Vl+^ 2 2 "l + ^ 2 3 J 



(C.6) 



(C.7) 



1 V^i+tf* 



In 



32 • 



d/3 ln(l + /3 2 ) 



y 00 d/3 ln( 

J-oo * 1 



and 



J-oo 7T y„ c 



12 



+ /3 2 



d/3i 



l+/3 2 2 l+/3 2 3 ; 2 l+/3 2 3 

(C.8) 



In 4. (C.9) 



In 4 



we obtain 



7T(l+/? 2 2 ) i-oo 27T 1 + (/3 12 -A) 2 1 + (^3 "A) 2 



^(P)^--ln4^--r'(P). 



(CIO) 



(C.11) 
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Expanding Ti in r' we obtain in the first order: 

x/ r dfh r dp 3 r°° dp, f i i \ 

^ ' loo 7T(1 + 7T(1 + loo 27T ^ + + 1 + fa) 

(cf ln(l + /3 2 2 ) + ln(l + (3D + ±Cf [ ln(l + /&) + ln(l + /&) ] ) 

(C.12) 

Using (C.9) and (CIO) we get: 

T^v) = 1 - \ In 4 Qcf } + Cf + Cf^ • r'(P) . (C.13) 

So, in the first order, 

T = T r {y) + Ti{p) = 1 - \ In 4 (2 Cf } + + ) • r\v) (C.14) 

References 

[1] S. Catani, L. Trentadue, G. Turnock and B.R. Webber, Nucl. Phys. B 407 (1993) 3. 

[2] S. Catani and B.R. Webber, Phys. Lett. B 427 (1998) 377 [hep-ph/9801350]. 

[3] S. Catani, G. Turnock and B.R. Webber, Phys. Lett. B 295 (1992) 269; 

Yu.L. Dokshitzer, A. Lucenti, G. Marchesini and G.P. Salam, J. High Energy Phys. 

01 (1998) Oil 

[hep-ph/9801324]. 

[4] V. Antonelli, M. Dasgupta and G.P. Salam, J. High Energy Phys. 02 (2000) 001 
[hep-ph/9912488]. 

[5] M. Beneke, Phys. Rep. 317 (1999) 1 [hep-ph/9807443]; 

B.R. Webber, Nucl. Phys. 71 (Proc. Suppl.) (1999) 66 [hep-ph/9712236]. 

[6] Yu.L. Dokshitzer, Perturbative QCD and Power Corrections, Invited talk at 11th 
Rencontres de Blois: Frontiers of Matter, Chateau de Blois, France, 28 Jun - 3 Jul 
1999 [hep-ph/9911299]. 

[7] B.I. Ermolayev and V.S. Fadin, Sov. Phys. JETP Lett. 33 (1981) 269; 
A.H. Mueller, Phys. Lett. B 104 (1981) 161. 

[8] Yu.L. Dokshitzer and S.I. Troyan, "Asymptotic freedom and Local Parton-Hadron Du- 
ality in QCD jet physics", Proceedings of 19th Leningrad Winter School, Leningrad 
1984, p. 144; for review see Yu.L. Dokshitzer, V.A. Khoze, A.H. Mueller and 
S.I. Troyan, "Basics of Perturbative QCD", ed. J. Tran Thanh Van, Editions 
Frontieres, Gif-sur-Yvette, 1991; 

A.H. Mueller, Nucl. Phys. B 213 (1983) 85, Erratum quoted ibid. B 241 (1984) 141; 
E.D. Malaza and B.R. Webber, Phys. Lett. B 149 (1984) 501. 



37 



[9] Yu.L. Dokshitzer, S.I. Troyan and V.A. Khoze, Sov. J. Nucl. Phys. 46 (1987) 712. 

[10] A. Banfi, G.P. Salam and G. Zanderighi, work in progress. 

[11] L3 Collaboration, B. Adeva et al., Z. Physik C 55 (1992) 39; 
D.P .Barber, et al., Phys. Lett. B 89 (1979) 139. 

[12] Yu.L. Dokshitzer, G. Marchesini and G.P. Salam, Eur. Phys. J. C 3 (1999) 1 
[hep-ph/9812487]. 

[13] A. Banfi, Yu.L. Dokshitzer, G. Marchesini and G. Zanderighi, under preparation. 
[14] S. Catani and M. Grazzini, [hep-ph/9908523] 

[15] A. Banfi, Yu.L. Dokshitzer, G. Marchesini and G. Zanderighi, under preparation. 

[16] Yu.L. Dokshitzer, A. Lucenti, G. Marchesini and G.P. Salam, Nucl. Phys. B 511 
(1998) 396 [hep-ph/9707532]. 

[17] S. Catani, G. Marchesini and B.R. Webber, Nucl. Phys. B 349 (1991) 635. 

[18] Yu.L. Dokshitzer, A. Lucenti, G. Marchesini and G.P. Salam, J. High Energy Phys. 
05 (1998) 003 
[hep-ph/9802381]; 

M. Dasgupta and B.R. Webber, J. High Energy Phys. 10 (1998) 001 [hep-ph/9809247]. 
[19] V.A. Khoze and W. Ochs, Int. J. Mod. Phys. A 12 (1997) 2949 [hep-ph/9701421]. 



38 



